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A framework is developed which quantifies the local exchange of energy and momentum between 
matter and the linearised gravitational field. We derive the unique gravitational energy-momentum 
tensor consistent with this description, and find that this tensor only exists in the harmonic gauge. 
Consequently, nearly all the gauge freedom of our framework is naturally and unavoidably removed. 
The gravitational energy-momentum tensor is then shown to have two exceptional properties: (a) 
it is gauge-invariant for gravitational plane-waves, (b) for arbitrary transverse-traceless fields, the 
energy-density is never negative, and the energy-flux is never spacelike. We analyse in detail the local 
gauge invariant energy-momentum transferred between the gravitational field and an infinitesimal 
point-source, and show that these invariants depend only on the transverse-traceless components of 
the field. As a result, we are led to a natural gauge-fixing program which at last renders the energy- 
momentum of the linear gravitational field completely unambiguous, and additionally ensures that 
gravitational energy is never negative nor flows faster than light. Finally, we calculate the energy- 
momentum content of gravitational plane-waves, the linearised Schwarzschild spacetime (extending 
to arbitrary static linear spacetimes) and the gravitational radiation outside two compact sources: 
a vibrating rod, and an equal-mass binary. 
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I. INTRODUCTION 

Half a century ago, a simple argument established that 
gravitational waves carry energy and can exchange this 
energy with matter. Often attributed to Feynman (cer- 
tainly popularised by Bondi [1]) the argument asked us 
to imagine a gravitational detector comprising a rigid rod 
along which two "sticky beads" are threaded. A passing 
gravitational wave then acts to alter the proper distance 
between the beads, and this motion, opposed by friction, 
heats the detector and thus mediates a transfer of energy 
from gravity to matter. Despite the simplicity of this 
idea, even after fifty years, it has not been possible to ex- 
plain where in spacetime this gravitational energy resides, 
and it is generally accepted that attempts to do so are 
"looking for the right answer to the wrong question" [2] . 

The elusiveness of the "right answer" , and the wrong- 
ness of the question, are very often identified as arising 
from gravity's gauge freedom, the consequence of which is 
a one-to-many mapping between physical spacetime and 
whatever localisation of gravitational energy-momentum 
might be proposed. Historically this issue was cast in 
terms of coordinate dependence, and the multitude of 
non-covariant objects that were constructed (first by Ein- 
stein [3], and most famously by Landau and Lifshitz [4]) 
were termed energy-momentum pseudotensors. However, 
a more recent formulation [5] has made it clear that the 
construction of a genuine tensor (defined on some back- 
ground spacetime) is not the central problem; rather, it is 
the tensor's dependence on the arbitrary diffcomorphism 
that maps physical spacetime to the background [6]. 
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Nevertheless, there is no reason a priori that gauge de- 
pendence should preclude the construction of a physically 
unambiguous tensor, provided we are prepared to remove 
the gauge freedom in some well-defined way. In cosmol- 
ogy this is frequently done by constructing new variables 
which are gauge invariant but equal to the relevant gauge- 
dependent fields (such as gravity or density fluctuations) 
in a particular gauge [7, 8]; however, it is just as ef- 
fective to provide a physically unambiguous method by 
which the gauge may be fixed, and to then insist that the 
gravitational field be evaluated in this gauge when locat- 
ing its energy and momentum. Unfortunately, no pre- 
vious approach has supplied instructions of this nature, 
and more importantly, neither the construction of these 
energy-momentum objects, nor their key properties, ap- 
pear to favour one gauge (or one set of gauge-invariants) 
over another; thus it appears impossible to justify any of 
these seemingly arbitrary choices as natural. 

Besides gauge dependence, there is also a great range 
of choice over which properties, physical or mathemat- 
ical, should define the gravitational energy-momentum 
tensor: should we be guided by a putative conservation 
law, or have in mind a particular role in the field equa- 
tions? For instance, it is always possible to locate the 
energy-momentum of matter by measuring the gravity 
it generates, so one might suggest that gravity's energy- 
momentum should be localised in a similar fashion, by ex- 
amining the interaction it has with itself. Following this 
idea to its conclusion, it has been shown [9-12] that gen- 
eral relativity may be constructed from an initially linear 
(spin-2) field theory that is then systematically coupled 
to its own (Hilbert) energy-momentum tensor. Sadly, 
this scheme leads us to identify the non-linear part of the 
Einstein tensor G a b — G^ as the gravitational energy- 
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momentum, so (a) the gauge problem remains, and (b) 
the result is additionally ambiguous, as different choices 
of "gravitational field" (g a b, g ab , \f~99 i eke-) mix the 
linear and non-linear terms in Gab- 
In spite of these various difficulties, one aspect of this 
enduring problem stands opposed to conventional wis- 
dom and motivates our present discussion: when gravity 
and matter interact, the exchange of energy is local! To 
see this we need look no further than the sticky bead 
detector: here, the energy exchange is certainly localised 
in so far as it takes place only within the confines of the 
detector. Furthermore, we can imagine a very small de- 
tector, much smaller than a wavelength of the incident 
gravitational radiation, and observe that at each instant 
a well-defined power is developed in the detector as heat; 
thus, at least in this case, the rate of energy exchange 
is associated with a particular point in spacctimc. One 
might hope, therefore, that consistency with this phe- 
nomenon would be enough to localise the energy and 
momentum of the gravitational field outside the detec- 
tor, or even when no detector is present. Moreover, even 
if a gravitational energy-momentum tensor could not be 
found, there would still be great value in constructing a 
framework for the description and analysis of local grav- 
itational energy-momentum exchange. The purpose of 
this article is to develop precisely this framework, and 
to examine the gravitational energy-momentum tensor 
it brings to light. In doing so we uncover a simple and 
unambiguous "right answer" through which the effects 
of gravitational energy-momentum may be usefully un- 
derstood. Conceivably, this was the "right question" to 
ask. 

For the sake of simplicity, we have restricted our 
present discussion to linearised general relativity on a fiat 
Minkowski background. It is only in this linear regime 
that the convenient fiction of a "gravitational field" prop- 
agating on a background spacetime can be taken se- 
riously, a construction which is essentially unavoidable 
when localising gravitational energy- momentum. 1 On a 
technical level, the restriction to the linear approximation 
limits the space of gauge transformations to a manage- 
able size, facilitating the analysis and eventual removal 
of our description's gauge dependence. Furthermore, our 
gravitational energy-momentum tensor will not be de- 



As long as there is some spacetime with everywhere vanishing 
gravitational energy-momentum, then this will naturally play the 
role of the background, and fluctuations away from this configu- 
ration will constitute the gravitational field. Although the most 
natural choice for this "ground-state" is flat spacetime, this does 
not necessarily preclude the extension of our formalism to less 
trivial backgrounds; however, we suspect there may be technical 
or conceptual problems with "ignoring" the energy-momentum of 
a nontrivial background. In particular, we anticipate issues anal- 
ogous to those of associating energy-momentum with a fluctua- 
tion in the electromagnetic field SF a i, when the background F a ^ is 
non-zero: the energy-momentum tensor T ~ F 2 + F8F + (SF) 2 , 
so the dominant contribution from the fluctuation will be linear 
in the field, rather than quadratic. 



rived from non-linear terms in the field equations, so we 
avoid any ambiguity arising from field redefinition. We 
shall not attempt to extend our results beyond the linear 
theory at this time. 2 

The structure of the paper will be as follows. We be- 
gin by building the foundations of our framework, de- 
riving a gravitational energy-momentum tensor (23) by 
demanding consistency with the energy and momentum 
exchanged with matter. As we will see, most of the ten- 
sor's gauge freedom is eliminated immediately as a nat- 
ural consequence of this derivation. Following this, we 
demonstrate two important additional properties of our 
tensor, further solidifying its interpretation as gravity's 
energy-momentum tensor. We then develop our frame- 
work more concretely by analysing the transfer of grav- 
itational energy-momentum onto an infinitesimal detec- 
tor; in the process of making this analysis gauge invari- 
ant, we will purge the last trace of gauge ambiguity from 
our energy- momentum tensor. Finally, we examine the 
gravitational energy-momentum in some specific exam- 
ples. Throughout, we work in units where c = 1, write 
hi = 8irG, and use the sign conventions of Wald [13]: the 
metric signature is (—,+,+,+), and the Riemann and 
Ricci tensors arc defined by [V c ,V d ]w Q = R a bcd v b , and 
Rab = R c acb- 

II. MOTIVATION AND DERIVATION 

The purpose of this section is to explain how, by con- 
sidering the energy-momentum transferred between mat- 
ter and gravity, we are led to a formula for the gravi- 
tational energy-momentum tensor. We begin by laying 
down some mathematical groundwork. 

A. Preliminaries 

As previously explained, this paper focuses exclusively 
on linear gravity: we only consider physical spacctimcs 
(Ai,g a b) in which the curvature R a bcd is everywhere 
small. As usual, this allows us to identify the physical 
spacetime with a flat background spacetime {M.,g a b), 
where R a bcd = 0, using a diffcomorphism (f> ; Jvi — > jCi- 
The "gravitational field" h ab is then defined on Ai by 

0*5ab = 9ab + h ab , (1) 

and we insist that <fi be chosen such that h ab is small 
everywhere, in order that the linearised Einstein field 



2 Of course, it may not be possible to extend the framework we 
develop here to the full non-linear theory, and we accept that 
localising gravitational energy-momentum in this regime (where 
the distinction between background and fluctuation is virtually 
meaningless) may be an inherently flawed idea. Of course, this 
does not alter the validity of our work in the linear case, where 
the "field theoretic" view is justified. 
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equations are a good approximation: 3 

G ab cd h cd = nf ab + 0(h 2 ). (2) 

In the above relation, T ab = 4>*T ab = 0(h) is the mat- 
ter energy-momentum tensor T ab mapped onto the back- 
ground, and 

G ab cd h cd = V c V {a h b) c - ±V 2 h ab - iV a V 6 /i 

+ hjab (V 2 h - V c V d h cd ) (3) 

is the linearised Einstein tensor G^ b . Our freedom of 
choice over cf> will of course give rise to the usual gauge 
transformation 6h ab = V( a ^)- 

On the background it will be useful to define four vec- 
tors 4 {e M a } obeying 

Vae£ = 0, (4) 

&[i &va = 1](U/: (5) 

which form the basis of a Lorentz coordinate system {x 11 } 
on Ai: e a = (d^) a . From this starting point, we shall 
define a corresponding set of vector fields {e^} in the 
physical spacetime, 

e/ = (^rV 1 ' (6) 

the behaviour of which will only be determined once we 
have fixed the gauge <f>, an issue to which we will return 
later. 



B. Energy-Momentum Currents 

Superficially, general relativity is a theory in which the 
energy and momentum of matter is always conserved: 

V a T ab = 0. (7) 

However, the sticky bead argument has already demon- 
strated that this is not the case; in reality, matter may 
gain (or lose) energy through interaction with the gravi- 
tational field. The reason for this apparent contradiction 
is as follows. In order to determine the energy of each 
part of the detector, one must first specify a timelike 
vector field e a (the "time direction" conjugate to the 
energy) with which to form an energy current-density 
ja — y» e i The incoming gravitational wave will then 
prevent e b from satisfying V a e & = 0, and we will find 
that V Q J a — T a b \7 a e b ^ 0. This inequality indicates a 



3 We use 0(h n ) as an abbreviation of 0({h ab ) n ); this should not 
be confused with the trace of the gravitational field h = h a t,g ab . 

4 We use Roman letters as abstract tensor indices [13, p. 437] and 
Greek letters as numerical indices running from to 3. Tensor 
indices of fields defined on the background are of course raised 
and lowered with g a i,. 



mismatch between the energy of the matter flowing into 
a given point, and the change in energy of the matter at 
that point; in other words, it represents the appearance 
of additional energy which was not already present in 
the matter - this is the energy absorbed from the gravi- 
tational wave! What is needed, therefore, is a framework 
which can account for this gained energy by identifying a 
corresponding loss in the energy of the gravitational field. 
We devote the rest of this section to the development of 
this idea, which will form the basis of our description of 
gravitational energy-momentum. 

Following the previous discussion, it should now be 
clear that we must define one energy current-density, and 
three momentum current-densities, by 

j/ = n e ;, (8) 

using the vectors {e M a } that get mapped to the Lorentz 
basis of the background. This is a generalisation of the 
practice of defining conserved currents by contracting T ab 
with a killing vector in a spacetime with a continuous 
symmetry. Here, however, the vector fields {e^"} only 
correspond to approximate symmetries (present because 
spacetime is nearly flat) and thus the currents will not 
be conserved. The real difficulty is choosing sensible be- 
haviour for {e^ a } that sufficiently captures the "paral- 
lelism" of killing vectors in the absence of any gravita- 
tional symmetry. Because e M a = (<fi~ 1 )*e fl a , this question 
has been recast as a choice of gauge, which we will ad- 
dress later. 

Having defined our energy-momentum currents (apart 
from specifying <p) we are now in a position to express the 
key idea of our approach. We seek a symmetric tensor 
field r ab , defined on the background, that is a quadratic 
function of the gravitational field h ab . We wish to be 
able to interpret r ab as the energy-momentum tensor of 
the gravitational field, and we shall achieve this by insist- 
ing that its non-conservation (in the background) exactly 
balances the non-conservation of the J^ a in the physical 
spacetime. Specifically, we wish to be able to define grav- 
itational energy-momentum current-densities j' a by 

J/ = r\e^\ (9) 

such that 

v4/ + 0*(v a j/)=o. (io) 

This equation captures the idea that energy-momentum 
is transferred between matter and the gravitational field. 
In particular, equation (10) indicates that knowing the 
behaviour of h ab at some point will be sufficient to de- 
termine the fields VaJ^ that express the local change 
in energy-momentum of the matter at the corresponding 
point in the physical spacetime. 

We proceed by calculating the two elements of (10). 
Because we are using a Lorentz basis in the background, 
V Q e^ = trivially gives 

Vai/ - eJ>V a r a b . (11) 
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The second term is a little less trivial; using (7), 

^(VaJ/) = 0*(nv aeAl b ), 

= ^2V*(V«e„ 6 ), 

= (f a b + 0(h 2 ))(V a e p b + e p c C b ac ), (12) 

where C a bc = ±(V b /i c a + V c h b a - \7 a h bc ) + 0{h 2 ) is 
the connection between the two derivative operators: 
(j)*(y a v h ) = W a (t)*v b +C b ac 4>*v c . Now, because V a e£ = 0, 
and T ab = T b(l7 we have 

0*(V„V) = i^&VflW + V ft c 6 - vXe) + 0(h 3 ) 
= ^ c T\V c h a b + 0(/i 3 ). (13) 

Finally, we use the field equations (2) to write 

0*(VaJ/) = ^e^{G ab cd h cd ){V q h ab ) + 0(h 3 ). (14) 

Inserting (11) and (14) into (10), and discarding the 
0(h 3 ) terms, we arrive at the defining relation of the 
gravitational energy-momentum tensor: 

KV a T aq = -\{V q h ab )G ab cd h cd . (15) 

The next step will be to use this equation to derive 
a formula for r ab in terms of h ab . In order to do so, 
however, we must make one additional demand: r ab will 
not depend on second derivatives of h ab , but will be a 
function of V c h ab and g ab only. The reason we must 
impose this condition is that equation (15) can only de- 
fine r ab up to the addition of "superpotential" terms, 
those fields whose divergence vanishes identically. Be- 
cause these terms are of the form S7 c S7 d H[ ac ][bd] (where 
H[ ac ][bd] = H [bd][ac] is some function of h ab ) they neces- 
sarily contain second derivatives; thus our restriction on 



T ab is sufficient to remove this ambiguity. At the mo- 
ment, it might be tempting to view this condition as a 
convenient way to tame the derivation, and keep in mind 
that we can always add in super-potentials later if we 
wish. However, in section III it will become clear that 
many of the interesting properties displayed by T a b will 
be unavoidably spoilt by the addition of such terms. For 
this reason we will not consider superpotcntials further 
here. 



C. Determining the Energy-Momentum Tensor 

In truth, it will not be possible to construct a sym- 
metric tensor r a b that satisfies (15) for all h a b; 5 to make 
progress we will need to impose some condition on V c h ab 
and specialise to this restricted set of gravitational fields. 
Although this forced restriction might appear to be a 
flaw in our formalism, as we shall soon see, it is actually 
a valuable asset. 

There are only three linear conditions we can place 
on Vc/ia& which neither introduce extra fields, nor break 
Lorentz invariance: (a) V c h ab = 0, (b) V a h = 0, or (c) 
V°/i a 6 = AVft/i, for some constant A. 6 Condition (a) is 
obviously far too restrictive: it does not allow us any 
gravitational field whatsoever. In contrast, condition (b) 
is not restrictive enough: there is no r ab that solves (15) 
for all gravitational fields with constant trace. 7 We must 
therefore focus on condition (c), which we repeat for later 
reference: 

V a h ab = \X7 b h. (16) 

Using this relation, it will be possible to replace any oc- 
currence of V a h ab with AV bh; hence the most general 
formula for a symmetric tensor r a b, & quadratic function 
of V c h a b, is as follows: 



KT pq = g pq (A V c h ab V c h ab + A 1 V a hV a h + A 2 V c h ab V b h ac ) + A 3 V p h ab V q h ab + A 4 V p hV q h + A 5 X7 a hV (p h q) a 

+ A 6 V a h b (p V q) h ab + A 7 V a h bp V a h q b + A s V b h ap V a h q b + A 9 V a hV a h pq , (17) 

where {A n } are arbitrary constants. We proceed by substituting this ansatz into (15) and solving for {A n }. First, let 
us calculate V p r pg by taking the divergence of (17); using (16) to convert every V a h a b to AV^ft., and collecting terms, 
we find that left-hand side of (15) amounts to 

KV p r pq = (2^o + A 3 )V q V c h ab V c h ab + (2Ai + Ai + ±AA 5 + XA 9 )V a hV a V q h + (2A 2 + \A 6 )V c h ab V q V b h ac 

+ A 3 V 2 h ab V q h ab + A±V 2 hV q h + {\A b + XA 7 + XA S + A 9 )V a V b hV a h bq + (|A 5 + \\A Q )V a V b hV q h ab 

+ \ArS a hV 2 h aq + (±A 6 + A 7 )V a V b h cq V a h bc + ±A 6 V a h bq V 2 h ab + A 8 V c h ab V a V b h cq . (18) 



5 We will shortly describe how to check this assertion, which is 
simply a property of (15) and independent of the requirement 
that T a b contain no second derivatives. 

6 The only other possibility, V a 7i a b = 0, can be achieved by taking 
(c) with A = 0. 

7 For the sake of brevity, we will not prove this assertion here. 



Instead we will attend to condition (c) and derive the formula 
for T a b that it admits. After we have done so, we invite the 
reader to perform a similar calculation under condition (b) and 
verify that no solution exists. 
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Meanwhile, (16) simplifies the right-hand side of (15): 

-±(V q h ab )G ab cd h cd = -\V q h ab ((A - ±)V a V 6 /i - ±V 2 /i ab + ig ab (l - A)V 2 /i) . 
Comparing (18) with (19), term by term, we conclude that the unique solution to (15) is 

A = -\, A 1 = ^ s , A 3 = |, Ai = ~l, A 2 = A 5 = Ae = A 7 = A 8 = Ag = 0, A = \. 



We have therefore determined the formula for our gravi- 
tational energy-momentum tensor, 

KT pq = ±\7 p h ab \7 q h ab - ±\7 p h\7 q h 

- \g Pq (y cK b V c h ab - \V a hV a h), (21) 

and the condition, 

V a h ab - \V b h = 0, (22) 

that we must place on the gravitational field. Finally, 
we introduce the abbreviation h ab ~ h ab — \g a bh for the 
trace-reversed gravitational field, and f a b = T a b — \sabT 
for the trace-reversed gravitational energy-momentum 
tensor, allowing us to compactly re-express our results: 

Kf pq = ±V p h ab V q h ab , (23) 

V a h ab = 0. (24) 

We are now in a position to justify our earlier claim, 
that the restriction on the gravitational field is not a hin- 
drance but a major advantage. The condition we have 
derived (24) is simply the defining equation of the har- 
monic gauged As this equation can always be satisfied 
by making a gauge transformation, it does not in any way 
limit the physical applicability of our approach! More- 
over, we have received a valuable gift: this condition 
has appeared as a natural consequence of the derivation, 
and forces upon us a very strong restriction for the dif- 
feomorphism </> that maps the physical spacetime onto 
the background. Specifically, </> _1 is required to map the 
Lorentz coordinates of the background onto harmonic co- 
ordinates in the physical spacetime. 9 We can therefore 
think of (24) as the condition that specifies the correct 
behaviour to demand of the basis {e^"} needed to de- 
fine sensible energy-momentum currents </ M a . We should 
stress, however, that while the harmonic gauge condi- 
tion has removed the vast majority of the gauge freedom, 
a small amount remains in the form of transformations 



This is also commonly referred to as de Donder gauge or Lorentz 
gauge. 

To see this, let {x^} be Lorentz coordinates on the background 
(V a V( ) a;' 1 = 0), and let {y^} be coordinates in physical space- 
time defined by J/ M (p) = x M (0(p)) for all p £ M. Then 
<£*(VV) = V 2 ^' - h ab V a ^J b x' 1 - CV a h ab - iVi,/i)V & ^ = 0. 



(19) 



(20) 

I 

Sh a h = V( a £h) which satisfy V 2 £ a = 0; we will return to 
this issue in section IV. 

This completes the derivation of r a b- We have found 
the unique symmetric tensor, a quadratic function of 
V c /i a b, that describes the transfer of energy and momen- 
tum between matter and the gravitational field accord- 
ing to (10). In doing so we have found that this solu- 
tion only exists if the harmonic condition V a /i ab = is 
obeyed, and this in turn has solidified the definition of 
energy-momentum currents J ' ° as the contraction of T ab 
with the basis vectors associated with harmonic coordi- 
nate systems of physical spacetime. In the next section 
we shall prove that r ab displays many other interesting 
properties very much in keeping with its interpretation as 
an energy-momentum tensor. In section IV we shall ex- 
amine energy-momentum exchange in detail, and address 
the last piece of gauge freedom. 



III. PROPERTIES 

Here we will demonstrate that, in two important spe- 
cial cases, T a h exhibits interesting mathematical proper- 
ties (beyond accounting for V a J^") that further promote 
its interpretation as the energy-momentum tensor of the 
linear gravitational field. 



A. Gauge Invariance of Plane- Waves 

Consider an arbitrary gravitational plane-wave: 

h ab = h ab {k^). (25) 

Here, k a is a constant vector, and {x 11 } are Lorentz coor- 
dinates on the background. The linear vacuum field equa- 
tion V 2 h a b = 0, and the harmonic condition \7 a h a b = 0, 
enforce 

k a k a = 0, k a h' ab = 0, (26) 

respectively, where the prime indicates differentiation 
with respect to the variable We wish to consider 

the most general gauge transformation Sh ab = V( a £(,) 
that maintains the plane- wave form of h ab - Clearly we 
require £ a = £ a (fc M :r M ), and thus 

Sh ab = k ia C b y (27) 
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Note that k a k a = now guarantees V 2 £ a = 0, ensuring 
that the harmonic condition (24) is not broken by the 
transformation. Let us now calculate the effect of this 
transformation on the gravitational energy-momentum 
tensor; working from (23), 

k5t„ = \V p 5h ab V q h ab + \V p 5h ab V q 5h ab 

+ \k p k q k ia ^r /b) -y ab k c o 

= l kp k q ((kr < !) 2 -(k c O 2 )=0. (28) 

Thus the energy-momentum of an arbitrary gravitational 
plane- wave is completely invariant under the gauge free- 
dom consistent with the harmonic condition and its 
plane-wave form. This is significant for a number of rea- 
sons. Firstly it reveals that, for the special case of plane- 
waves, we need not concern ourselves with the gauge free- 
dom that remains after enforcing the harmonic condition: 
the requirement that the gauge be chosen such that the 
plane-wave form of the field be manifest is sufficient to 
unambiguously define the energy-momentum tensor r ab 
from the physical spacetimc (M.,g). Thus, even if one 
does not accept our method for resolving the last of the 
gauge ambiguity (to be presented in section IV) it is still 
possible to stop at this point and agree that a well-defined 
energy-momentum tensor for gravitational plane-waves 
has been found. Secondly, this particular gauge invari- 
ance will prove useful when we wish to produce a global 
picture of the motion of energy-momentum: if the source 
region of a gravitational wave is very far from the detec- 
tion region, we may use a different gauge in each and yet 
still produce a consistent picture of energy-momentum 
transfer - the field in intermediate region will approxi- 
mate a plane-wave, and thus r a b in this region will agree 
with both end-point gauges. 10 There is also a third signif- 
icance to this result, but this will only become apparent 
once we have demonstrated the second important prop- 
erty of gravitational energy-momentum: positivity. 



B. Positivity 

This section concerns the energy-momentum of 
transverse-traceless (tt) gravitational fields, those for 
which h = 0, \7 a h ab = 0, and u a h a b = 0, for some 
constant timclike vector field u a defined on the back- 
ground. 11 For now let us simply suppose that these con- 
ditions apply to h a b and derive the consequences for r a b- 
We shall justify our interest in this specialisation, and 
offer an interpretation of u a , in section IV. Presently, let 
it suffice to say that because these conditions may always 



10 This idea is explained fully in section IV D. 

11 Clearly, the harmonic condition is satisfied as a result of these 
requirements. 



be imposed (at least locally) by a gauge transformation 
in regions where T ab = 0, the results of this section will be 
generally applicable to vacuum regions, but it will not be 
necessary to demand that T ab = globally. We present 
our result as the following theorem. 

Theorem. //, at some point p € M., the gravitational 
field h ab obeys the transverse-traceless conditions 

V a h ab = 0, h = 0, u a h ab = 0, (29) 

for some timelike vector u a , then r a b satisfies the follow- 
ing inequalities 

v a T ab v b > 0, (30) 
v a T ac r c b v b < 0, (31) 

at p, for any timelike vector v a . 

Proof. Without loss of generality we can set u a u a = — I 
and v a v a = — 1. Now, introduce two Lorentzian coordi- 
nate systems at p, the first {x°,x 1 } (i = 1, 2, 3) such that 
u° = 1, u* = 0, and the second {x ' ,x 1 '} (i' = 1,2,3) 
such that v° = 1, v l = 0. The transverse-traceless con- 
ditions (29) reduce (23) to 

KT pq = -\g pq V c h ab V c h ab + \W p h ab W q h ab , (32) 

and set h^i = hoo = 0. Using the primed basis to express 
the tensor indices of g ab and V a , the unprimed basis for 
h ab , and writing h%j = do'hij, we find that 

KV a T ab V b = KTo'O' 

= \m,h l3 f-{k J ) 2 ) + \{h l3 Y 

= !((^<%) 2 + (M 2 )>0, (33) 

because (di>h i:i ) 2 = J2i,j,i>=idi'hijdi'hij and (%) 2 = 
Si j=i do'hijdo'hij are sums of squares. Similarly, 

(Wr ab ) 2 = I6 K 2 (-(t , ,) 2 + (r 0V ) 2 ) 
= -l{{ dl ,h lJ f + {h lJ ff 
+ {h i , J d il hij){h kl di>h k i) 

= -Kidi.hijf - (kjff - (M 2 (c^M 2 

+ {hijdi'h^^had^hki) 

= -\m>h tJ ) 2 - (^) 2 ) 2 

- \0iijdehki - hud^h^) 2 < 0. (34) 

□ 

The inequalities we have just deduced are the gravi- 
tational version of the Dominant Energy Condition: the 
first indicates that r ab only ever defines positive energy- 
densities; the second indicates that the flux of this energy 
can never be spacelike. Succinctly, they tell us that grav- 
itational energy is positive and never flows faster than 
light. As the Dominant Energy Condition has always 
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referred to matter, we will avoid confusion if we resist 
subsuming (30) and (31) under this name; instead, when 
the gravitational energy-momentum tensor obeys these 
inequalities (for all timelike v a ) we shall simply say that 
it is positive, and write r a f, > shorthand. 

That T a i, is positive for all transverse-traceless h a b is 
one of the major advantages our approach has over pre- 
vious descriptions of gravitational energy-momentum [3- 
5, 14, 15]. Provided we work with a transverse-traceless 
field (with respect to some timelike vector u a ), which is 
always possible locally in a vacuum, t & will always make 
good physical sense in that it will obey its own version 
of the Dominant Energy Condition. To some extent, this 
result supplies its own justification for choosing the tt- 
gauge whenever possible; however, we will see in the next 
section that these conditions arise naturally by consid- 
ering the gauge invariant transfer of energy-momentum 
onto point-sources. Furthermore, the significance of u a 
(in terms of energy-momentum transfer) will also be ex- 
plored through these arguments. 

Before we move on, however, we take this opportu- 
nity to present an important corollary of the plane- wave 
gauge-invariance of section III A. It is well known that 
there always exists exactly one gauge transformation of 
the form (27) that takes an arbitrary plane-wave (25) 
to one obeying the TT-conditions [16]. Hence we can 
transform any gravitational plane-wave into transverse- 
traceless gauge without altering the energy-momentum 
tensor, at which point the positivity theorem ensures 
that T a b > 0. Thus, all gravitational plane- waves have 
positive energy-momentum tensors, even if they are not 
transverse-traceless. 



IV. INTERACTIONS 

In this section we apply our formula for the gravita- 
tional energy-momentum tensor to the interaction be- 
tween gravity and an idealised matter distribution that 
we shall refer to as a point-source. The reader will be fa- 
miliar with the compact source, an isolated body confined 
to a compact spatial region of radius d much smaller than 
the wavelength A of the gravitational radiation it emits; 
point-sources are the limit of such systems as d — > 0, 
entirely analogous to the infinitesimal dipoles of electro- 
magnetism. 12 Not only will this provide a useful example 
of the practical application of our approach, these con- 
siderations will finally allow us to rid ourselves of the last 
trace of gauge-dependence in our description. 

From now on we will work almost exclusively in the 
flat background spacetime; as such it will generally be 
convenient to represent all tensors in some Lorentzian 
coordinate system {x^}, and to drop the "caron" mark 



from T^y. Thus, our formula for the gravitational energy- 
momentum tensor is written as 

Kf M „ = jdf.ha/jdyh^ , (35) 

the harmonic condition becomes 

0"V - 0, (36) 

and the linearised field equations (2) are 

d 2 V = -2kT„ v . (37) 

Also, it will be useful to separate {x^} into a time coor- 
dinate t = x°, and spatial coordinates x = (x , x 2 , a; 3 ), 
define a radial coordinate r = and to use lower-case 
Roman indices {i, j,k . . .) to indicate spatial components. 
Typically, the coordinates will be implicitly chosen to co- 
incide with the rest-frame of the system under consider- 
ation. 



A. Pulses and Point-Sources 

The core of our analysis will be to examine the most lo- 
calised gravitational interaction possible: an infinitesimal 
point-source (at x = 0) met by an instantaneous pulse 
plane-wave (propagating along the x 1 direction, arriving 
at x = at t = to). 

As a result of calculations in appendix A, we know 
that a point-source has the following energy-momentum 
tensor: 

Too = MS(x) + ^IijdidjSix), 

T oi = |(4 + Jtf)0j*(2), (38) 

Tij = \h 3 5{x). 

Here M and J,j = Jr^-i are constants representing, 
respectively, the mass and angular momentum of the 
source, and ly = Iuj\(t) its (time dependent) quadrupole 
moment. 13 We do not intend to use this point-source as 
an actual source of gravitational radiation, but rather as 
a probe of the energy-momentum of the incident pulse. 
To this end, we are interested in the limit M, Iij , Jy — > 0, 
allowing us to neglect the self-interaction of the source. 
As this procedure is entirely analogous to using a "test- 
particle" to probe the geometry of spacetime, we shall 
refer to the point-source as a test-source in this limit. 

The gravitational field will consist of two parts: h^ u = 
^source + ^wavc ThfJ first; due to the tes t-source, is given 

in appendix A by (All) and satisfies the inhomogeneous 
field equations 

^"source = _ 2kT ^. (39) 



We derive the energy-momentum tensor and gravitational field 
of the point-source in appendix A. 



The reader should refer to appendix A for definitions of these 
quantities in terms of the infinitesimal limit of the compact 
source. 
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The latter is the incident pulse plane- wave, 



^-™ vc — A fiU H(k a x a — to), 



(40) 



where H is the Heaviside step function, k^ = (1, —1, 0, 0) 
is a null vector in the x 1 direction, and A^ v is a con- 
stant tensor satisfying k^A^ v = 0, as demanded by the 
harmonic condition. Obviously, h^ ve satisfies the homo- 
geneous field equations: <9 2 /i™ ve = 0. 

Let us now compute 9 m t m „, which, via (10), quanti- 
fies the exchange of energy-momentum between the test- 
source and the gravitational wave. Starting from (35), 
we have 

K d^ v = \d 2 h a pd v h a V 

= -\KT a pd v {A^H{k a x a -t )) 

+ o((^T cc ) 2 ), (41) 

and we neglect terms of order (h s °^ rcc ) 2 compared to 
those of order h™ ve h s ™ IC0 in the limit M./y, Jy -> 0. 14 
Using H' = 5, the Dirac delta function, we arrive at 

9%, = ~kkj{k a x a - t o )T a0 A a ? 
= -\k v 8{k tT x' T -t ) 

x (lijAijS(x) - 2(iij + J i:j )dj5(x)A i0 



(2MS(x) + I lJ d l d J 5{x))A 00 



(42) 



This is the equation we sought. It determines the energy 
and momentum collected by our probe due to the incident 
pulse, and locates this transfer in spacetime. The key 
problem is that above relation is not, as it stands, gauge 
invariant; we address this issue the next section. 



B. Gauge Invariance and Microaveraging 

The incident wave possesses gauge freedom that nei- 
ther breaks the harmonic condition nor spoils its pulse 
plane-wave form: 



5h™ ve = d( M £„) , 



^ = ^A(fc Q x Q -t ), 



(43) 



where is any constant vector, and A' = H. The 
effect of this transformation is to alter A^ by 5A^ V = 
k( p E v }, and although the transverse components do not 
change (SA22 = SA 23 = 8A 33 = 0) the right-hand side 



14 There is a slight technical issue here. From (All), we can see 
that, as r -> 0, h B °^ TCC -4 00; thus h s °^ TCC inevitably becomes 
larger than /i™J ve at small enough distances. Strictly speaking, 
then, one should use a finite-size source (of radius d, say) when 
one takes M,I,J -s- and neglects 0((/i=°„ urce ) 2 ). As we can 
choose d to be as small as we like, however, we can always replace 
the finite source with an equivalent point-source after this limit 
has been taken. 



of (42) is clearly not invariant. The beauty of working 
with an instantaneous interaction, however, is that we 
can average over the (infinitesimal) interaction region 

lim B e (t ), 

where B € (to) = {(*, ■ I* - *o| < e, W | < e}, (44) 

without sacrificing the localised description of 9 m t m „ . Let 
us call this operation a microaverage (at x = 0, t = to) 
and denote it by (. . .) to : 



{f) to = 5{x)5{t-t ) 



lim 

r->0 



fd 4 x. 



Be(io) 



(45) 



For the interaction we are analysing, the integral 
Ib (t ) d^T^v^x captures the key physical content of 
d^T^i,. To elaborate: the divergence theorem equates 
this integral with J QIS , 4 . r^d^S^ which measures the 
mismatch between the net flux of gravitational energy- 
momentum entering through a spherical surface S e = 
{x : \x\ < e} barely larger the source, and the gravi- 
tational energy-momentum contained within S e that is 
gained between the times t = to — e and t = to + e. By 
the defining property (15) of t^ Vi this mismatch in grav- 
itational energy-momentum precisely accounts for the 
energy-momentum absorbed by the source, which is what 
we wanted to know. The only information we have lost 
in taking the microaverage is the knowledge of precisely 
where within the test-source the energy-momentum is be- 
ing absorbed. As we have let the size of this probe shrink 
to zero, however, this is of little concern. 

The computational advantage of the microaverage is 
that the integration in (45) allows us to transfer deriva- 
tives off the delta- functions in (42); for example, 



5(k a x a — to)iijdjS(x)Aiod 4 x 



(46) 



Be(to) 

djS(t — x 1 — to)iij5(x)Aiod i x 
S(t — x 1 — t )i i i8(x)A i Qd A x 

B t {t ) 

= - I 8{t-x l - t )/;i<5(x)A j0 d 4 x 

JB t {t ) 
= —Iil(to)Aio- 

Applying this technique to the whole of (42) yields 
x (iijAij + 2i a A i0 + IuAqo + 2MA 00 

Finally, we unpack k^A^ = 0, 

Aqo + An + 2A i = 0, A 22 + A 33 = 0, 
A02 + A 12 = 0, A 03 + A 13 = 0, 



(47) 



(48) 
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and substitute these into (47). The result is 

(& l T lu ,) to = -$k v 6(g)6(t-t ) 

x (l x A x +I+A + + MA oo y (49) 

where we have written the transverse components of the 
wave as A x = A23 and A + = (A22 — A^)/2, and 
extended this notation to Iij. We are almost done: 
SA X = 5A+ = under the gauge transformation (43), 
so the first two terms in (49) are manifestly gauge invari- 
ant; however, the term proportional to MAqq is not. 

Various arguments can be made to show that this 
"monopole term" is physically irrelevant to the energy- 
momentum transfer we are considering. At the simplest 
level, the fact that we are free to set Aqo to any value (in- 
cluding zero) through gauge transformation (leaving A x 
and A + untouched) indicates that the monopole term 
can have no bearing on the energy-momentum of the 
physical system under scrutiny. Furthermore, if we con- 
sider a wave for which A^ = k(aE u \, then it is clear 
that, while such a pulse is gauge-equivalent to flat space- 
time (Afj, v — 0) it would nonetheless register a transfer 
of energy-momentum if the monopole term were to be 
believed. 

The physical irrelevance of the monopole term should 
come as no great surprise, as there can be no way to ex- 
tract energy-momentum from a gravitational wave using 
a monopole alone (i.e. a test-source with = Jjj = 0): 
an observer sitting on an isolated point mass could per- 
form no local test to distinguish whether a gravitational 
wave had even passed, and in particular, must be unable 
to extract any energy. 

In fact, all that the monopole term is responding to 
is a change in normalisation of the time coordinate in 
the physical spacetime: </>*(egeo a ) = — 1 + h o- Naively, 
we might expect this factor to be significant as it repre- 
sents the Newtonian potential at the test-source. How- 
ever, this is not a local effect. The only way an observer 
on the test-source could be aware of such a shift is by 
comparison with some standard clocks at spatial infin- 
ity. The pulse plane-wave prevents this idea from be- 
ing well-defined, however, as it divides spatial infinity 
into two regions: a; 1 < t — to, where wave has already 
been received, and x 1 > t — to, where it has not. For- 
tunately, a gauge can always be chosen that does not 
suffer from this inconsistency; setting Aqo = is the only 
way to ensure that the standard clocks at infinity all run 
at the same rate (relative to our coordinate t) and this 
inevitably removes all trace of the monopole term from 
the interaction. Thus the insistence that the clocks at 
infinity agree with each other amounts to a prescription 
that removes the gauge-dependence of our microaveraged 
energy-momentum transfer. 15 We can implement this 
procedure mathematically (without fixing the gauge, or 



15 The reader should not be under the impression that the monopole 



setting M = 0, which is physically untenable) by acting 
on {d^T^ v ) to with the operator (1 — M8m)' 

(d»T ia ,)g = (l-Md M )(ff>T^) t0 

= -±kj(x)5(t-t ) (l x A x +I + A + ^) . (50) 

We shall call this the monopole-free microaverage. This 
is a local, completely gauge invariant description of the 
energy-momentum transferred onto test-sources by pulse 
plane-waves. Furthermore, the right-hand side of (50) 
has an obvious physical interpretation: the coupling be- 
tween Iij and Aij can be understood, roughly speaking, 
as the product of a force (responsible for accelerating 
the constituents of the quadrupole moment) and a dis- 
tance (actually an expansion/contraction of spacetime) 
and thus represents the work done on the test-source. 
For example, consider a test-source composed of two bod- 
ies of mass m separated by a light elastic rod of length 
2d aligned with the x 2 -axis; provided the amplitude of 
the motion of the masses is much smaller than d, then 
I22 — 4dma, where a is the (outward) acceleration of 
each mass. Due to the gravitational wave, the proper 
distance of each mass from the centre of the rod in- 
creases by ^22^/2; thus, counting the motion of both 
ends of the rod, the total work done on the source, by 
the wave, is -2(ma)(A 22 d/2) = -I 2 2^22/4 = -I+A+/2, 
which agrees precisely with (50). Thus we see that the 
monopole-free microaverage corresponds to the familiar 
physical quantities that we would intuitively use to define 
the energy and momentum of the test-source. 

In the next section we will generalise the monopole- 
free microaverage to arbitrary gravitational fields, and 
uncover a substantial mathematical shortcut that will 
greatly simplify this procedure. 

C. Arbitrary Gravitational Fields 

Clearly, pulse waves are a special case, and one might 
expect that for an arbitrary (harmonic gauge) plane- wave 

h™ ve = B^(k a x a ), WB UV - 0, (51) 

we would need to perform finite averages, rather than 
microaverages, to remove the gauge dependence of our 
description; thus, we would be forced to sacrifice our lo- 
calised picture of energy-momentum transfer. However, 



term is universally insignificant. Thus far we have argued its 
irrelevance only for pulse plane-wave, and as we shall see at the 
beginning of the next section, this idea follows by linearity to 
general gravitational waves. However, should the gravitational 
field have a time-independent part, then it is possible for this to 
couple to the monopole in a physically meaningful way. This is 
due to the particularly limited gauge freedom available to hoo 
when the field is time-invariant. We will return to this issue in 
Section IV E. 
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7 wave 



provided B^ u (t) — > as t — > — oo, we can always write 
B liv {t {) )5(k a x a - t )dt 

B^{t )H{k a x a -t )dt 
[B^(t )H(k a x a -t )]tZ 

B^(t )H(k a x a -t )dt , (52) 



and perform the monopole-free microaverage on each 
component of this sum: 

(3%,[/C/T + B a p(to)H(k„ar - t )])f o dt . 

(53) 

The result of this process is 



(0%„)f = -\K5{x) (l x B x + I+B 



(54) 



which renders the interaction completely gauge invariant, 
and does not sacrifice the local character of our descrip- 
tion of energy- momentum transfer. 16 

Although the operation of splitting the wave into a se- 
ries of pulses and performing a monopole-free microaver- 
age on each pulse may seem too complicated to be useful, 
the same result can be achieved by a simple alternative 
method: transform h™^ vc to transverse-traceless gauge, 
where the vector u' 1 referred to by the TT-conditions (29) 
corresponds to the rest-frame of the test-source. Then, 
when we calculate d^r^, we will automatically recover 
the monopole-free microaveraged result. To demonstrate 
this, we recalculate d^T^ vl generalising (42) for use with 
arbitrary plane-waves (51), 



+ (2MS(x) + Iijdidjdix^Boo) , (55) 



and substitute the TT-conditions B 0u = B = (which, 
along with k^B^ v = 0, set B\ v = and B 2 2 = — -B33): 



d»T™ = d' J 'T^[h sonice + (/i wave ) TT ] 
= —-^k v 8{x)IijBij 
= -\K5{x) (l x B x +I+B A 



Hence, 



(56) 



(57) 



Just as we write Iij for Iij(t), we have, in (54) and elsewhere, 
left the argument of B MJ/ (i) implicit. 



Furthermore, this equation is not only applicable to inci- 
dent plane- waves. Because both sides are linear in h^ vc , 
equations (57) must also hold when h™^ vc is any sum 
of plane-waves, propagating in arbitrary directions. Lo- 
cally, we can always express h™^ vc as a sum of plane- waves 
(and some time- independent part, which we will ignore 
until section IV E) so, quite generally, we have 



= -\8(x)iijd v K 



TT 
>J ' 



(58) 



where hj^, is the incident gravitational field in transverse- 
traceless gauge. This equation provides an easy method 
for calculating the energy-momentum transferred onto 
the microaveraged test-source due to the presence of ar- 
bitrary incident gravitational radiation. Moreover, we 
see that (ignoring the time-independent field) our gauge 
invariant probe only exchanges energy-momentum with 
the transverse-traceless field; the other components of the 
field do not play a role in this process. We explore the 
wider significance of this result in the next section. 



D. Energy-Momentum and Transverse- Traceless 
Gauge 

Let us now take a step back from the fine details of the 
test-source interaction and assess the general picture that 
is unfolding. As we first saw in section II C, is not in 
general invariant under the gauge freedom that remains 
after the harmonic condition has been enforced. As a 
partial remedy of this, the monopole-free microaveraged 
test-source emerged as a local, gauge invariant probe of 
gravitational energy-momentum exchange. It has now 
come to light that only the transverse-traceless field takes 
part in this process. From this standpoint, a method sug- 
gests itself which will remove the remaining ambiguity of 
t m „ in a natural fashion: simply transform the incident 
gravitational field to transverse-traceless gauge! Conse- 
quently, only the degrees of freedom relevant to gauge 
invariant energy-momentum exchange will contribute to 
the gravitational energy-momentum tensor. We shall 
codify this idea as a gauge-fixing program defined in 
terms of two "frames" , one associated with gravitational 
detectors, the other with astrophysical sources. In what 
follows, the gauge-fixing only refers to the dynamical part 
of the gravitational field; as we explain in section IV E, 
the time- independent part of the gravitational field is es- 
sentially gauge invariant and so does not need to be fixed 
in any way. 

Detector- frame. Consider a gravitational detector D 
in a region Vd which contains no matter besides the de- 
tector. We shall suppose that the incident field (due to 
sources outside Vd) is much larger than the field due 
to the detector itself; in other words, we model D as 
a test-source. The detector-frame is then obtained by 
transforming the incident field to TT-gauge, taking u M 
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to be the four- velocity of the detector. 17 As a result, 
the energy-momentum transferred onto D will be exactly 
equal to the gauge invariant quantities defined by the 
monopole-free microaverage. What is more, we can imag- 
ine adding hypothetical test-sources (co-moving with D) 
anywhere within Vd in order to "measure" the gravi- 
tational energy-momentum there; because the field has 
been prepared in this gauge, the result will agree with 
the gauge invariants we have already defined. In this 
way, the detector-frame defines t^ v through the gauge 
invariant energy-momentum that would be absorbed by 
furnishing Vd with an array of infinitesimal probes mov- 
ing at the same velocity as the actual detector. 

Source-frame. Now consider a compact source S in 
a region Vs which, as above, contains no other matter. 
In contrast to the detector, we shall assume any inci- 
dent field can be neglected in comparison to the outgoing 
field due to S. The source-frame is obtained by trans- 
forming the outgoing field to TT-gaugc, taking to be 
the four-velocity of the source. This gauge transforma- 
tion can only be achieved by breaking the harmonic con- 
dition at S (see Appendix B for details) so it will not 
be possible to use t^ v to describe the energy-momentum 
lost by the source; however, this self-interaction is ill- 
defined for a point-like system anyway, and we would 
have to resolve the source into component parts before 
such a question could be answered. Outside the source, 
h^ v will remain harmonic, so t m „ will still represent the 
energy-momentum that could be absorbed by a hypo- 
thetical test-source (with the velocity of S) were we to 
insert one. Much like the detector-frame, we can think 
of this prescription as measuring t^ v by filling Vs with 
infinitesimal probes that are co-moving with the source. 

The gauge-fixing program is simple: if one wishes to 
describe the energy-momentum of the gravitational field 
as it would be measured by some detector D, then adopt 
the source-frame near distant astrophysical sources, and 
the detector-frame of D everywhere else. This allows the 
energy-momentum in the vicinity of a source to be unam- 
biguously determined by the source alone, whilst simulta- 
neously adapting the gravitational field for a description 
of energy-momentum absorption by the detector. Re- 
markably, despite using different TT-gauges in the vari- 
ous regions, this program still produces a self-consistent 
picture of the propagation of energy-momentum from 
the sources to the detector. This is because, many 
wavelengths from an isolated source, the gravitational 
field approximates a plane-wave, and so (as we saw 
in section III A) the gravitational energy- momentum is 



17 Note that it is not the total field, but just the incident field, 
which is made transverse-traceless. We cannot alter the gauge 
of the field generated by the detector because it is impossible to 
produce outgoing spherical waves in the gauge field £ M without 
breaking the harmonic condition at D. This cannot be allowed to 
happen if Tuv is to account for the energy-momentum exchanged 
with the detector. 



:s 2 . 



Source-frame 

Plane-waves 

Detector-frame 



FIG. 1. Schematic showing plane-wave regions sewing var- 
ious source-frames to the detector-frame of D. Si and S2 
represent compact sources many wavelengths apart; if they 
have different velocities, then each will determine a separate 
source-frame. In contrast, 53 and S4 are only separated by a 
small number of wavelengths; as there is no plane- wave region 
dividing the sources, they must share their source-frame. It is 
natural to base this joint source-frame on the velocity of the 
centre of mass of the multi-component system. 



gauge invariant there. Thus, in this regime, the source- 
frame energy-momentum and the detector-frame energy- 
momentum are equal, regardless of the relative velocity 
between the source and the detector. The plane-wave re- 
gions will therefore "sew together" the different frames 
and produce a globally consistent description of energy- 
momentum flowing from source to detector. When the 
sources are not isolated, but arc separated from each 
other by only a small number of wavelengths, no plane- 
wave region will exist between the sources. In this case, 
the sources must be treated as one extended source, with 
a joint source-frame that identifies with the four- 
velocity of the centre of mass of the many-body system. 
We illustrate the gauge-fixing program schematically in 
figure 1. 

Up to this point, we have justified our insistence 
on TT-gauge based purely on considerations of energy- 
momentum exchange with matter. Of course, there is 
another exceptional property of t^ v , derived in section 
IIIB, that also holds under these conditions: it is always 
positive. This is a peculiar and surprising result. It is 
something of a small miracle that transverse-tracelessness 
guarantees not only agreement with the monopolc-frec 
microaverage, but also ensures that represents posi- 
tive energy-density, and causal energy flux. Under these 
conditions it will always be possible to make physical 
sense of the gravitational energy-momentum tensor: we 
will never have to interpret (or explain) negative or su- 
perluminal energy. 

We feel that the dual significance of transverse- 
traceless gauge leaves little doubt that this is the correct 
procedure by which to remove the final trace of ambigu- 
ity in the definition of Tav . In section V we will apply this 
program to a small number of examples, including two 
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specific compact sources: a vibrating rod, and an equal- 
mass binary. First, however, we must address a technical 
issue regarding the time-independent part of the gravita- 
tional field. 



E. Time-Independent Fields 

Unlike the dynamical part of the gravitational field, 
the gauge of the time-independent mode (i.e. the time- 
averaged field) is completely fixed by insisting that (a) 
the harmonic condition d^huv = holds everywhere, 
(b) h^y — > as r — > oo, and (c) gauge transformations 
dhfj,i> = are finite everywhere and bounded at in- 

finity. To see this, suppose that our time-independent 
field h^ v (x) obeys the harmonic condition and vanishes 
at spatial infinity. Then the transformed field h' = 
h^v + d(u£v) will only satisfy the harmonic condition if 
<9 2 £ = 0, and will only be time- independent if £ p = 0. 
Thus df(d^^)) — 0, the only bounded solutions of which 
are constants (by Liouville's theorem). Hence we are 
forced to take h!^ v = if the new field is to also van- 
ish at spatial infinity, and we thus conclude the gauge is 
unique. 18 

This result reveals that we are not required to perform 
any form of microaverage to remove the gauge depen- 
dence of the energy-momentum transfer associated with 
the time-independent mode of the gravitational field: this 
mode is already gauge invariant. In truth, this is a rather 
convenient situation. We could not microaverage a time- 
independent field even if we needed to, due to the caveat 
B^uit) — > as t — > — oo, encountered when deriving (53). 

As there is no gauge freedom in the time-independent 
part of the gravitational field, we cannot expect this 
mode to have h afi = or h = in general. This leaves 
open the possibility (at least in principal) that close to 
sources, where the time-independent mode can become 
comparable in amplitude to the dynamical field, the pos- 
itivity of may be compromised. However, in sec- 
tion VC we will see that, even very close to (but not in- 
side) a compact source, the time-independent field obeys 
h 00 ^> h 0i » hij. 19 It is easy to show that such a field 
will not upset the positivity of r Ml/ . Neglecting the small 



This argument relies on our insistence that the harmonic con- 
dition be valid everywhere. In section VC we will make use 
of the following mathematical trick: by relaxing the harmonic 
condition at the source itself, we will be able to combine many 
local TT-gauges to form a gauge in which the dynamical part of 
the gravitational field (outside the source) is transverse-traceless 
for all time. As we will see, however, it is impossible to ap- 
ply this procedure to a time-independent field. Thus there is 
nothing to gain from weakening the harmonic condition on the 
time-independent mode, and it is therefore kept unbroken. 
These order-of-magnitude inequalities arc not limited to the com- 
pact source. The time-independent mode of the gravitational 
field is always generated by the time-averaged energy-momentum 
tensor of matter (T M „ ) , and it is to be expected that this field will 
be dominated by the slow (i.e. non-relativistic) motion of matter, 



quantities, the trace-reversed gravitational field will take 
the form 

h IMV = K^ v -±§u IJ ,'u v , ( 59 ) 

where K^ v is the dynamical field in transverse-traceless 
gauge, and $ = — hm/A is the the Newtonian poten- 
tial, the only non-negligible contribution from the time- 
independent mode. The trace-reversed gravitational 
energy-momentum tensor therefore takes the form 

iKf^ = d ll h (X f S d v h a0 = d^hjjdvhjj + 8d^$, (60) 

which ensures that the positivity proof of section IIIB 
can proceed almost exactly as before, with $ effectively 
behaving as an additional component of hJJ . Thus, 
even though it is not transverse-traceless, the time- 
independent mode does not give rise to any negative or 
superluminal energy. 

V. APPLICATIONS 

This section is devoted to calculating the energy- 
momentum content of the gravitational field in a small 
number of examples, following the gauge- fixing program 
of section IV D. 



A. Plane-Waves 

Although we have already studied gravitational plane- 
waves in a variety of contexts, we have yet to evaluate 
the energy-momentum they carry. This calculation will 
serve as a simple first example, and will illustrate the use 
of the detector-frame. 

We begin with an arbitrary (harmonic gauge) plane- 
wave, 

h^ = h^(k a x a ), fc"fc„ = 0, k»h'^ = 0, (61) 

and substitute this field into equation (35): 

KT^v = ^k^k v h a jjh a ^ . (62) 

As we ascertained in section III A, the energy-momentum 
of plane- waves is gauge invariant. Consequently, we can 
simplify (62) by evaluating in TT-gauge, thereby re- 
moving all components except for h + and h x - 

KV = fak v ((ti + ) 2 + (ti x ) 2 ). (63) 

In fact, because 5h + = 5h x =0 under any gauge trans- 
formation that keeps a plane-wave, the right-hand 



so that (Too) > (T 0i ) > (Ty>. Hence h 00 > h 0i > h tJ will hold 
quite generally: whenever (T^) is dominated by non-relativistic 
motion. 
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side of this equation is gauge invariant also. Hence, equa- 
tion (63) must hold in any gauge, and all other terms 
on the right-hand side of (62) must cancel in general. 20 
Using this formula for r Mt/ , every future-directed time- 
like unit- vector defines a gravitational energy current- 
density, 



v%„ = v^k^K{{h') 2 + (h' x ) 2 )/2 K , 



(64) 



which is clearly future-directed and null; unsurprisingly, 
the energy of a gravitational plane-wave is positive and 
flows at the speed of light in the direction of propagation. 

So far, the gauge invariance of t m „ has made gauge- 
fixing unnecessary. The insistence that we evaluate the 
/i M „ in the detector-frame only becomes important when 
there are multiple plane-waves propagating in different 
directions. Suppose, for example, that there are two 
plane- waves: 



(65) 



where k\ and k 1 ^ are non-parallel null vectors. As is 
quadratic in h^ v , the energy- momentum of the total field 
takes the form 



where r* and rj^, are the individual energy-momentum 
tensors of hr„ v and /ij 1 , respectively. Now, any gauge 
transformation that preserves the form (65) of the grav- 
itational field can be thought of as a pair of gauge- 
transformations that act on h}- and h l ^ v separately, pre- 
serving their plane- wave forms; thus t^ v and rS, must be 
invariant under gauge-transformations of this type. How- 
ever, the "cross-terms" in (66) are gauge-dependent, as 
the (gauge-dependent) longitudinal components of ft* 
will be transverse to ft?„, and vice versa. This gauge am- 
biguity is removed, however, by the presence of a physical 
detector: once we demand that the energy-momentum 
exchanged with this detector is to equal the monopolc- 
free microaverage, we fix the gauge completely. This is 
the detector-frame: is transverse-traceless, with w M 
identified as the four-velocity of the detector. In this 
sense, the gauge-fixing program is the procedure that 
enables us to "add together" the energy-momentum ten- 
sors of gravitational plane- waves (which, individually, are 
gauge invariant) to form the energy-momentum tensor of 
the total field. 

For the sake of the concreteness, let us set fc^ = 
(1,1,0,0), k 1 ^ = (1,0,1,0), and u* 1 = (1,0,0,0). Then, 
once we have transformed ft Mt , to TT-gauge, the energy- 
momentum tensor becomes 



ii/ 



kt^v — kt^j, + kt 1 ^ v — ^{2k l ^k 1 ^ + 7j flL ,)h 1 + h + 



This can be verified by taking k^ 
k^h^y = in the same form as (48). 



(1,-1,0,0) and using 



where h\ 



n 22 



-h\ 3 , and ft+ 



h 11 

"■33 



-ft 11 



Due 



to the positivity theorem of section IIIB, we already 
know this tensor describes a positive energy-density, and 
a causal energy flux. As a particular example of this, it is 
easy to calculate the energy-density associated with 



«r 00 = I((^) 2 + (^) 2 + (^') 2 + (O 2 ) 



±h v h lv 

4' t +' t + ! 



(67) 



(ft+') 2 > h l lh l l'/2, we can confirm that 



and, as (ft+) t yu>+ 
this energy- density can never be negative 



B. Linearised Schwarzschild Spacetime 

The Schwarzschild spacetime is the vacuum solution 
to the Einstein field equations outside any uncharged 
spherical non-rotating body of mass M. At distances 
much greater than kM, where the linear approximation 
is valid, the gravitational field must therefore correspond 
to that of the compact source with Iij = Jij = 0: 



too 



kM 
2irr ' 



(68) 



( 66 ) and h 0i 



0. Obviously, this is an example of a 



gravitational field that is entirely time-independent; thus, 
as explained in section IV E, there will be no possibility of 
transforming to transverse-traceless gauge, nor any need 
to do so. 21 The formula (35) for the gravitational energy- 
momentum tensor yields 



/ M 



\87rr 2 



(69) 



where i M is the radial unit vector. 22 It is easy to confirm 
that this energy-momentum tensor is everywhere posi- 
tive. Any timelike unit vector defines a positive grav- 
itational energy-density, 



(2(^) 2 + 1 )> ' (70) 



and an energy current-density J" 
nowhere spacclikc: 



J V J V 



-K 2 



M 



< 0. 



which is 



(71) 



It is also worth comparing equation (69) with the elec- 
tromagnetic energy-momentum outside a point-charge: 



21 Of course, the linearised Schwarzschild spacetime can be rep- 



resented in other gauges, but no others obey d^h^ 



and 



h^v = everywhere, and are well-behaved at infinity. 

Xi/r from ap- 



22 This is a trivial extension of the notation x 
pendix A; we simply define xq = 0. 
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(flyii/ + 2u ll u u — 2x^x u )/r A . While both tensors that ho^ = and h = continues to be true for t £ (to 



diminish in proportion to 1/r 4 , they define very differ- 
ent stress profiles at each point. The gravitational field 



has 



-Tee 



tqo > and thus describes ra- 



dial compression, tangential tension, and negative pres- 
sure; while the electromagnetic field has — T rr = Tee = 
Ta = Too > and thus describes radial tension, tangen- 
tial compression, and positive pressure. The physical sig- 
nificance of this difference is far from obvious, but may 
relate to the like-attracts-like character of gravity: the 
negative gravitational pressure mediating the attraction 
of other masses, while positive electromagnetic pressure 
causes the repulsion of like-charges. In addition, it may 
be possible to understand the radial gravitational com- 
pression (and tangential tension) in terms of some "elas- 
tic" analogy for spacetime, as the Schwarzschild geome- 
try "squeezes in" extra radial distance (between spheres 
of given area) in comparison to fiat space. However, the 
theoretical value of such an analogy is unclear, and we 
do not intend to develop it any further here. 

Although we have focused here on the linearised 
Schwarzschild spacetime as a particular example of a 
time-independent field, we note in passing that it is 
easy to evaluate the gravitational energy-momentum ten- 
sor associated with any static configuration of matter 
T^i, = u^u v p(x): equation (35) simplifies to 



(72) 



where the Newtonian potential $ = — hoo/i is deter- 
mined by solving Poisson's equation = Kp/2. Equa- 
tion (72) reveals that the energy-momentum of the New- 
tonian potential is exactly that of a massless Klein- 
Gordon scalar field. 



C. Gravitational Field of a Compact Source 

We shall now calculate the energy-momentum content 
of the gravitational field (All) generated by a compact 
source. 23 The first step will be to enter the source-frame: 
we must transform the dynamical part of the outgoing 
field into TT-gauge, with u M identified as the four- velocity 
of the source. We can always make this transformation 
locally by choosing the gauge fields £ M such that h 0fl = 
hou = and h = h = at some time t = to; then 



r, to + r). This method is problematic in that it is based 
around an arbitrary special time irji arL d that transverse- 
tracelessness always breaks down within a time At = 2r; 
these issues prevent us from forming a global picture of 
the energy-momentum outside the source. 

As we show in appendix B, these problems can be 
completely avoided if we weaken the harmonic condition 
slightly, so that d^h^ = is only enforced outside the 
source. This trick allows us to find a gauge in which the 
dynamical field is transverse-traceless everywhere outside 
the source, for all t, and does not require us to choose 
a special time to- We can think of this gauge as a way 
of joining up the many possible local gauges (defined us- 
ing the aforementioned method) in a mutually consistent 
fashion. 25 The process of transforming the gravitational 
field of the compact-source (All) can be found in ap- 
pendix B, here we simply display the result: 



h 00 = {2M + (I^didj) 



hoi — Jijdj-^— 



H 

47rr 



2tt 



hiSij + IStkSji - 4Jk(i5j)ij dkdi 



LO 2 



uj 2 ( ISij - 2I l3 



8nr 



where we have introduced the notation 



(Iij) ee lim / I^t) 



(73) 



-dt, (74) 



for the time-average of the quadrupolc moment, and 



[i l3 {t)-(h ))&t, 



(75) 



for the Fourier transform of its dynamical part. Notice 
that the terms proportional to M, Jy, and (Iij) consti- 
tute the time-independent mode of the field, and have 
therefore not been transformed. At this point we can 
confirm the assertion of section IV E, that the time in- 
d' 2 ^ = (which preserves the harmonic condition) and dependent field satisfies h 00 > h 0i > hy. To do so 
the field equations d 2 h^ u = (outside the source) ensure 



This calculation should not be confused with the analysis per- 
formed in section IV, where a test-source (essentially a compact 
source in the limit d, M, Jij , Iij — > 0) interacted with an incident 
field, which presumably had been generated by another source, 
very far way. Here the compact source will represent an astro- 
physical source (with finite d, M, Jij and Iij) and by adopting 
the source-frame we will compute the energy-momentum of the 
outgoing field as it would be measured by microaveraged detec- 
tors co-moving with the source. 



See [13, chap. 4.4b] for details. 

Presumably, there is some topological obstruction which prevents 
us from joining these local gauges without violating the harmonic 
condition at the source. However, provided we do not intend 
to calculate the energy-momentum transferred between matter 
and gravity at the source, this is not an issue. Even if we were 
careful to keep the gravitational field harmonic at the source, 
more work would be needed to perform such a calculation, as 
this self-interaction only becomes well-defined by breaking down 
the compact source into component parts. 
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we note that, firstly, there is no time- independent term 
in hij, and secondly, seeing as the radius of the source 
d > J/M, and that we are outside the source (which is to 
say, r 3> d, the regime of validity of (73)) then we must 
have J/r < M. 

Having rendered the dynamical field transverse- 
traceless outside the source, all that remains is to sub- 
stitute (73) into (35) to calculate t^. As was shown 
in the process of deriving (60), the energy- momentum 
of the time-independent field adds linearly (i.e. without 
cross-terms) to that of the dynamical field. Given that 
we have already investigated the part due to the time- 
independent field in section VB, it is generally more 
interesting to discard this term, and focus on the ad- 
ditional energy-momentum due to the dynamical field. 
In figure 2 we show the results of a computation of this 
additional gravitational energy-density Too outside two 
monochromatic compact sources: a vibrating rod, and 
an equal-mass binary. It goes without saying that the 
energy-density is everywhere positive, and that the en- 
ergy current-density is nowhere spacelike. 



VI. CONCLUSION 

It is natural to suspect that wherever matter gains 
energy under the influence of gravity, a corresponding 
loss in the energy of the gravitational field must have 
occurred. By constructing a framework to quantify this 
idea, we have succeeded in localising the energy and mo- 
mentum of the linear gravitational field, and have shown 
this energy to be positive and to not flow faster than 
light. 

The core result of our investigation is the formula 
(23) for the gravitational energy-momentum tensor, the 
unique symmetric tensor, quadratic in V c /i a 6, which ac- 
counts for the energy-momentum lost or gained by mat- 
ter through its interaction with gravity (10). Crucially, a 
tensor satisfying these conditions only exists in the har- 
monic gauge (24) and thus, as a matter of necessity rather 
than choice, our framework discards nearly all its gauge 
freedom. A small set of viable gauge transformations still 
remain, however, and although these do not alter the 
energy-momentum of gravitational plane- waves (§111 A) 
this invariance does not extend to arbitrary gravitational 
fields. 

In response to this ambiguity, the monopole-free mi- 
croaveragc was developed (§IVB); this constitutes a 
local and fully gauge-invariant description of energy- 
momentum transfer, and agrees with the intuitive no- 
tion that the "work done" on a gravitational detector is 
the product of the force (proper acceleration) and the 
proper distance through which the force is applied. Of 
the incident field, only the transverse-traceless part con- 
tributes to the microaveraged exchange (57), and thus a 
natural gauge-fixing program is motivated, based around 
transverse-traceless gauge (§IVD). The effect of this 
program is to prepare the field so that no microaver- 



age is needed, and furthermore, to ensure that energy- 
momentum is only assigned to those components of the 
field whose energy-momentum can be measured by a 
microaveraged detector. Because the positivity prop- 
erty (§111 B) holds true wherever the field is transverse- 
traceless, the gauge-fixing procedure also ensures that 
(for the dynamical field at least) gravitational energy- 
density is positive, and gravitational energy flux is time- 
like or null. No longer burdened by gauge ambiguity, 
the gravitational energy-momentum tensor can be eval- 
uated without difficulty: the energy-momentum of grav- 
itational plane- waves (63), the linearised Schwarzschild 
spacetime (69), and the gravitational radiation outside 
compact sources (Fig. 2) have been provided as specific 
examples. 

With regards to further investigation, there are two 
obvious directions in which our framework might be ex- 
tended: beyond the linear approximation, and beyond 
the flat background. 26 However, it is currently unknown 
whether such extensions are possible, or even conceptu- 
ally sound. On a more practical level, one could apply 
our formalism to the energetics of actual gravitational de- 
tectors, rather than the idealised test-sources that have 
so far dominated our discussion. In doing so, the frame- 
work developed here may benefit the design and analysis 
of future gravitational-wave experiments. 
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Appendix A: Sources 

The aim of this Appendix is to derive the formula 
for T^u that defines a gravitational point-source (essen- 
tially an infinitesimal gravitational quadrupole) and the 
field h^y that it generates. The derivation comprises two 
parts: first, a calculation of the field due to a compact 
source; second, a calculation of the field due to a candi- 
date T^u that vanishes everywhere but at x = 0. As the 
field from the first calculation matches that of the sec- 
ond (within the region of validity of the compact source 
approximation) we will be able to conclude that our can- 



To extend T a (, beyond the linear regime, one would hope to con- 
struct a tensor t ab , defined on the physical spacetime M, such 
that !p*tab — T ab + 0(h 3 ). Clearly, it will only be possible to 
make this identification if T a i is gauge invariant to second or- 
der, as S((f>*t a i,) ~ 0(h 2 )dt; under a change of gauge, whereas 
Sr a i, ~ 0(/i)<9£ unless it is invariant. Thus, only once r a j, has 
been brought into the detector-frame, or the source-frame, can 
we proceed. Consequently, wc should expect that t ab will not 
only depend on the physical metric g a i,, but also on the four- 
velocity of the relevant detector or source. 
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FIG. 2. Plots of the energy-density of the dynamical gravitational field outside two monochromatic compact sources: a vibrating 
rod, and an equal-mass binary. Only half a period is shown, as t m „ oscillates with twice the frequency of the source. Although 
the rod and the binary are much smaller than one wavelength, they have been magnified to illustrate the phase of their 
motion. The propagation of gravitational energy is more easily appreciated in the animated versions of these plots, available 
at www.mrao.cam.ac.uk/~lmb62/animations . 



didatc Tpy is indeed the energy-momentum tensor we 
sought, that of an infinitesimal compact source. 



1. The Compact Source 

A compact source is an isolated gravitational body con- 
fined to a compact spatial region T> of radius d much 
smaller than the wavelength A of the gravitational radia- 
tion it emits. Although calculations of the field h^(x, t) 
outside a compact source are available in many stan- 
dard references, we present our own here for two rea- 
sons. Firstly, textbook treatments commonly conflate 
the slow-motion approximation (d -C A) with the far- 
field approximation = r A). Here we shall assume 
only that the source is very small (d <C r, A) but not any- 



thing about the ratio of A to r. 27 Secondly, the standard 
approaches frequently omit a full calculation of /loo and 
hoi- Presumably, these components are ignored because 
they do not appear to contribute to the gravitational field 
in transverse-traceless gauge; however, they must be in- 
cluded if is to satisfy the harmonic condition. 

The retarded solution to the linearised field equations 
(37) is given by 



K f T^ v {x , t \x 



\x — x' 



(Al) 



We shall proceed by expanding the right-hand side of this 
equation to second order in the small quantities d/X and 



27 As we are working with in the confines of linearised gravity, we 
should also insist that d 2> 2kM, the Schwarzschild radius of the 
source. However, this will have little bearing on our calculation. 
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d/r, so that we have an integral of energy-momentum 
tensors T^ v = T„ v (x' ,t — r) evaluated at the same time 
t' = t — r. Using 

I* " *\ = r (l - + |f - + 0((d/r)»)) , 

(A2) 



equation (Al) expands to 
h, 



K 

lixr 



dV 



T)xv I 1 



r z T, 



X ■ X 

x ■ x 1 

{ x-x'f 
2r 4 



1\2 



2r 2 

'12 



2r 2 



3(f-f) 

2r 4 
3(x-£') 2 

2^ 



0((d/r) 3 ) 



(A3) 



Although we have not written their arguments, it should 
be understood that the T^ v terms in the integral are eval- 
uated at (x*,t — r), while is evaluated at (x,t). 

In order to relate this integral to the basic physical 
properties of the source, we define its mass, momentum, 
and dipole moment by 



M= I T 00 dV, 



Pi = 

Xi = 



T 0l d 3 x', 



TooXjd x , 



(A4) 



respectively. Notice that, because the source is entirely 
contained within T> (so = on the boundary dT>) the 
conservation equation d^T^ v = (the linearised version 
of (7)) leads to the following relations: 



Xi 



P, 



daToox'A 3 x' 



(c^oK-dV 



- / r j0 (^4)dV = Pi, 



l ,A\ \r' = 0. 



(A5) 
(A6) 



Thus Xi = 0, and we are free to fix Xi = P. ; = by 
our choice of coordinate system. Note also that M = 
follows by an identical argument. Next we define the 
quadrupolc moment 



rp ' ' j3 / 

TooXiXA x , 



V 



and then derive 



la = -2 / To^dV, 

T> 

I T„dV, 



(A7) 

(A8) 
(A9) 



in a similar fashion. Finally we define the angular mo- 
mentum of the source 



(A10) 



and note that conservation sets J,j = 0. 

Before substituting these definitions and results into 
(A3), note that equations (A8) and (A9) indicate that 
/ ToidV - I/d - Md/X and jT^x' ~ I ~ Md 2 /X 2 ; 
hence the integrals of Tqj and Tij already have (respec- 
tively) one and two extra factors of (d/X) than the in- 
tegrals of Too. Thus, to second order, hij will include 
contributions from only the zeroth order quantities mul- 
tiplying in (A3), and hoi will include only first and 
zeroth order quantities multiplying To;. The final result, 
accurate to second order in the small quantities {d/X) 
and (d/r), is therefore 



SlijXiXj 




(All) 



where Xi = Xi/r is the radial unit vector, and all the hj 
terms are evaluated at the retarded time t' = t — r. Note 
that, while the fields /loo and hoi are often omitted from 
standard calculations, even in the far-field limit (r — > oo), 
they still contain terms of equal size to hij] these are 
necessary for consistency with the harmonic condition. 

We have successfully derived the form of the gravita- 
tional field outside a compact source. However, because 
(All) was constructed under the approximation scheme 
d <C r, we can only trust these equations at distances 
much larger than the size of the source. However, we can 
still ask the following question: what source would pro- 
duce a field such that (All) was valid for all r, no matter 
how small? This is the point-source we have been inter- 
ested in: the limit of the compact source as d — > 0. In the 
next section we present a candidate for the point-source, 
calculate its gravitational field, and show that this agrees 
with (All) for all r. 



The Point-Source 



Consider the following energy-momentum tensor for 



matter: 



Too = MS(x) + ilijdidj5(x), 
To* = \{iij + Jij)dj5{x), 



(A12) 
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where M, Jij — are constants. 1^ = I(ij){t) is inde- 
pendent of x, and overdots indicate differentiation with 
respect to t. It is easy to check that this distribution 
obeys d^T^ = 0. 

We wish to solve the linearised field equations 

d 2 V = -2kT^, (A13) 

looking for the retarded solution. Recalling that 

d 2 (f(t - r)/r) = -47r6(x)f(t), (A14) 

for any twice differentiable function f(t), we see that we 
can replace / — > nlij / 4ir to generate the result 



nljj (t - r) 
4nr 



(A15) 



Also, from equation (A14), we have 

&(dj(f(t - r)/r)) = -4^f(t)dj5(0). (A16) 
Thus, setting / — > n{jij + Jij) /4m gives 
' Iij(t - r) + \ 



h 0i = -^d 3 



-in 



By the same method, 



(A17) 



Iij(t- 



i - kM -+- K a a 

2irr 47r \ r 
k / 

4tt \ r r 2 



(A18) 



Therefore the source (A12) generates a gravitational field 
identical to that of the compact source (All), except that 
these equations are now valid for all x (except, possibly, 
x = 0) not just r 3> d. The energy-momentum tensor 
(A12) is the point-source we required and (All) the field 
it generates; the correspondence with the compact source 
allows us to validate the interpretation of M as the mass, 
Iij the quadrupole moment, and J^ the angular momen- 
tum of the source. 



Appendix B: Persistent Transverse- Traceless Gauge 

Here we describe a method by which the dynamical 
part of the gravitational field outside a compact source 
(centred at x = 0) may be transformed to a gauge which 
remains transverse-traceless for all time, everywhere out- 
side the source. This will be achieved by relaxing the har- 
monic condition slightly, so that d^h^ = only holds 
outside the source. 



First, a point of notation. The gauge transformation 
described in this section is only applicable to the dynam- 
ical part of the gravitational field h^J 1 = — (h^), 
where (...) signifies a time average. Rather than crowd 
the notation, it will be convenient to assume that (h^) = 
0, and use /i M „ to stand for h^J 1 . For the compact source, 
this amounts to setting M = J^ = (ly) = in (All). 
At the end of the calculation we will reinsert these time- 
independent terms to the transformed field without al- 
teration. 

The general procedure is as follows. To begin, take the 
Fourier transform of the dynamical part of the gravita- 
tional field: 



h^{uj,x) = 



t h^ iv {t,x)&t. 



(Bl) 



The Fourier transform renders the field equations as 

(uj 2 + df)h^(u;,x) = 0, (B2) 

everywhere outside the source, i.e. for i / 0. The har- 
monic condition becomes 



iu}h n + dihin = 0, 



(B3) 



and writing £ M (w, x) for the Fourier transform of ^(t, x), 
the general gauge transformation 8h^ v — d<^ u ) takes the 
form 



Shoo = M; , 

6h 0i = ^iujii + 5 Oifo i 

5hij = 

To achieve transverse-tracelessness we set 

Co = iu)~ l hoo, 

ii = 2iuj~ 1 h M - uj~ 2 dih 00 - 



(B4) 



(B5) 



From the field equations (B2) it is clear that this gauge 
transformation obeys (w 2 +df)£ fl — for x ^ 0, and thus 
the harmonic condition is preserved outside the source. 
It is also easy to check that (B5) fixes Shoo = —hoo and 
Shoi — ~hoh an d hence ensures that the transformed 
field h' = h^ v + 8h^ v has h' 0fl = everywhere. Further- 
more, 

6h = —5h Q + Shu 

= h Q0 + di{2iuj~ 1 h Ql - uj~ 2 dih 00 ) 

= -uj~ 2 d 2 h 00 - ha, (B6) 

where, in the last step, we have used the fi = compo- 
nent of (B3). Thus, for f^0, where we may use (B2), 
we have 



5h = hoo — hu 



(B7) 
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Standard Method 



Persistent TT-gaugc 




= 0, 






h = 0. 


Is 




to 


Si 












V 





- 0, h „ ^ 0, h ^ 0. 



V, ± 0, /i „ ^ 0, h ^ 0. 



FIG. 3. Comparison of the standard method for achieving 
transverse-traceless gauge in the vicinity of a source [13, chap. 
4.4b], and the "persistent" method described here. In the 
two diagrams, S represents an arbitrary source (a region with 
Tf_i V ^ 0) moving relative to u M . The hypersurface t — to used 
to define the gauge in the standard method is also shown, but 
plays no role in our method. 



so that h! = outside of the source. In summary, the 
transformed field is 



/>'; 



(B8) 



with all other components zero, and h! = 0, d^h' = 
everywhere outside the source. 28 

We are now in a position to apply this procedure to the 
gravitational field of the compact source (All). Before 
doing so, however, it is worth mentioning that the tech- 
nique just described is not limited to compact sources. In 
generalising, the only adjustment needed is that (B2) will 
only hold at x such that T^(t, x) = for all t. Figure 3 
illustrates the difference between this technique and the 
standard method mentioned in section VC. 

Continuing with the compact source, we write the dy- 



namical part of (All) as 



hoi = d {ni t:j {t - r)/4irr), 
hij = Klij(t - r)/47rr, 



(B9) 



and take the Fourier transform: 



/i 00 = Iijdidjine %ur /Anr), 
/4irr, 



(B10) 



where Uj is the Fourier transform of the dynamical part 
of the quadrupole moment. Substituting this into (B8) 
yields 



Kj = 
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IkiSij + IS lk Sji - 4/ fe (j<5j)H d k di 



US 2 



U 2 ( I8« - 21; 



8irr 



(Bll) 



Finally we recall that h' = h'^ u (for x ^ 0) and reinsert 
the time- independent mode 



(h 00 ) = (2M + (I^didj) —, 



Airr 



(hoi 



(B12) 



(ha) = 0, 
to confirm equation (73). 



It should now be clear why this method cannot be applied to 
the time-independent mode of the field: ill-defined contributions 
proportional to S(uj)/u) or 8(lj)/uj 2 would appear in the integral 
on the right-hand side of (B8). Even without a delta- function 
at w = 0, this integral is not unambiguous until we explain how 
to deform the contour to avoid the poles there. We suggest the 
contour should dodge into the lower half of the complex plane, as 
this ensures that (ti ) is dependent only on h^v (t2 ) for < t\ , 
which is to say, the transformed field does not depend on future 
values of the untransformed field. Using this "causal" contour, 
we can substitute (Bl) into (B8) and perform the ui integral, ar- 
riving at = hij{t,x) + f*^ dt'((t - t')didjhoo(t' , x) - 
2d(ihj) Q (t' , x)). In general, this formula is less useful than (B8), 
however it does reveal the asymptotic conditions that the dy- 
namical field must obey for this gauge-transformation to be well- 
defined: as t — > —oo, the non-oscillatory modes of didjhoo and 
9(jhjjQ must vanish faster than t~ 2 and t _1 respectively. 
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